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We propose a model of suppression of 7r°'s based on two different effects: at low px we 
take into account the shadowing corrections, which are non-linear and essential for the 
description of the inclusive spectra, while at large pt the suppression is produced through 
the interaction of the large px pion with the dense medium created in the collision. The 
main features of the data on AuAu and dAu collisions at RHIC energies are reproduced 
both at mid and at forward rapidities. 
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1. Introduction 

One of the most intriguing results of the heavy ion program at RHIC is the fact 
that the nuclear modification factor in central AuAu collisions, RAA{b,y-,PT) = 
dyd^pT / n (b) JydPpT ' ^ s smauer than 1 for all pt- That means that the yield 
of particles produced in AA collisions at mid-rapidities and large px increases with 
centrality much slower than the number of binary collisions n(b). This phenomenon 
is particularly interesting since it is not observed in dAu collisions at RHIC at 
mid-rapidities. 

In order to explain these data, one should take into account two kinds of effects. 
On one side, we have the effects that are important at low px- In most of the mod- 
els of multiparticle production, two contributions to the multiplicity are considered: 
one proportional to the number of participant nucleons, N part , and a second one 
proportional to the number of inelastic nucleon-nucleon collisions, N^, t , which is 
the dominant contribution at asymptotic energies. In order to get the right multi- 
plicities at RHIC it is necessary to lower this contribution. Different mechanisms 
have been proposed, like saturation in high density QCD through the Colour Glass 
Condensate ^ or clustering of strings in the string models. Another possibility is 
the shadowing through pomeron interaction in the Dual Parton Model 1^1 All these 
mechanisms have in common the modification of a multiple scattering pattern -in 
the target rest frame- or gluon interaction -in a fast moving frame-. At very high 
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energies, the shadowing corrections -which are non- linear- lead, as s —* oo, to a sat- 
uration of the distributions of partons in the colliding nuclei. These corrections are 
important for the description of inclusive spectra, for particles with px ^< px >■ 
But with increasing px the shadowing corrections decrease and the scaling with 
n(b) is predicted in perturbative QCD (Raa — > 1)- 

So in order to explain the data at large pt, it is neccessary to consider that for 
particles with large momentum transfer, there are in general final state interactions 
-jet quenching or jet absorptionl^L. These interactions lead to an energy loss of the 
large pt parton (particle) in the dense medium produced in the collision. Hadrons 
loose a finite fraction of their longitudinal momentum due to multiple scattering. 
From this point of view, it is natural to expect that a particle or a parton scattered 
at some non zero angle will also loose a fraction of its transverse momentum due 
to final state interactions with a medium of partonic or prehadronic nature. The 
scattered particle does not disappear as a result of the interaction but its px is 
shifted to smaller values. Due to the steepness of the px distribution, the effect of 
suppression at large px may be quite large. Moreover, in this case there is also a 
gain of the yield at small px due to particles produced at larger px - which have 
experienced a px shift due to the interaction with the medium. This suppression 
vanishes at low px- when the px of the particle is close to < px >, its px can 
either increase or decrease as a result of the interaction, i.e. in average the px shift 
tends to zero. Our aim here is to describe the suppression of the yield of pions in a 
framework based on final state interactions and taking into account shadowing and 
Cronin effect. 



2. The model 

Our approach contains dynamical, non linear shadowing. It is determined in terms 
of diffractive cross sections and it applies to soft and hard processes. The reduction 
of multiplicity from shadowing corrections can be expressed by 

Jd 2 sf A (s)f B (b- S ) T A {b) 

RMb) = Tab(s) ' Mb) = T+AFWW) W 

where the function F(s) represents the integral of the triple Pomeron cross section 
over the single Pomeron one: 

ry m ** 1 d 2 a ppp 
F(s)=4tt dy | f =o = C[exp(y max ) - exp(y min )]. (2) 

J ymin a \s) dydt 

From this formula one can observe that for a particle produced at y = 0, with 
Vmax = h^ n ( s / m x) an d Vmin = ln(i? J 4mjv/ V3) , the effect of shadowing decreases 
at large px due to the increase of mx in y m ax- 

In other words, due to coherence conditions, shadowing effects for partons take 
place at very small x, x <C x cr = 1/jtinRa where rriN is the nucleon mass and Ra 
is the radius of the nucleus. Partons which produce a state with transverse mass 
mx and a given value of Feynman xp, have x = x± = x \ + 4m^/s ± if). 
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Our shadowing applies to soft and hard processes. Nevertheless, for large px these 
effects are important only at very high energies, when x ~ ^ satisfies the above 
condition. At fixed initial energy (s) the condition for existence of shadowing will 
not be satisfied at large transverse momenta: In the central rapidity region (y* = 0) 
at RHIC and for of jets (particles) above 5(2) GeV/c the condition for shadowing 
is not satisfied and these effects are absent. 

So for the large px partons it is neccessary to include the final state interactions. 
The interaction of a large px particle with the medium is described by the gain and 
loss differential equations which govern final state interactions: 

-n- = -o- Ps [Ph{x,Pt) - Ph(x,p t + dp T )\ (3) 

where ps and p B correspond to the space-time density of the medium and of large 
Pt 7r°'s, <j is the interaction cross-section averaged over momenta and d 4 x represent 
the cylindrical space-time variables: longitudinal proper time r, space-time rapidity 
y, and transverse coordinate s. 

Assuming a decrease of densities with proper time 1 /r -isoentropic longitudinal 
expansion, transverse expansion is neglected- the above equation transforms into: 

t dN -< b ^y^ = ^ N(bi s , y ) [ Njr0 (6, 8 , y, PT ) - Ko (b,s,y, PT + 8p T )\ (4) 

where N(b,s,y) = dN/dy d 2 s(y,b) is the density of the medium per unit rapidity 
and per unit of transverse area at fixed impact parameter, integrated over px, and 
N^o(b, s, y, P x) is the same quantity for 7r°'s at fixed px- 

If we integrate from initial time tq to freeze-out time Tf and taking into ac- 
count the inverse proportionality between proper time and densities, t//to = 
N(b, s, y)/N pp (y), where N(b, s, y) is the density produced in the primary collisions 
in DPM and N pp (y) is the density per unit rapidity for minimum bias pp collisions 
at y/s = 200 GcV= 2.24 fm~ 2 , we obtain the suppression factor S v o(b, y,pr) of the 
yield of 7r°'s at given px and at each impact parameter, due to its interaction with 
the dense medium. We get: 

c ^ / d 2 s a AB (b) n(b,s)S 7r o{b,s,y 1 p T ) 

J d 2 s &ab [b) n(b,s) 
where the survival probability is given by: 

Nn»{b,s,y,p T + Sp T ) 



S v o(b,s,y,p T ) = expi -a 



1 



N w o(b,s,y,p T ) 

Here a AB (b) = {1 - exp[-a pp AB T AB (b)}}, T AB {b) = J d 2 sT A (s)T B (b - s), T A {b) 
are the profile functions and n(6, s) = AB a pp T A (s) T B (b — s)/a AB (b). 

One can estimate which amount of the effect takes place in the partonic phase 
and which one happens in the hadronic phase. We can divide our suppression factor 
in two parts: 

Partonic, from initial density N(b, s, y) = d ^ R [i V ~ to d ^i V ~ ^p-j or equiv- 
alently from tq = 1 frn to r p = 3.36 fm, and hadronic, from partonic density 
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d^lpL ~ Jgo. t0 N ( y ) = = 2.24 fm- 2 , or equivalent^ from r„ = 3.36 fm 

to r/ = 5 — 7 fm. We find that 75% of the effect takes place in the partonic phase 
while only 25% of the effect takes place in the hadronic phase. 



3. Numerical results 

In order to perform numerical calculations, we need a -our free parameter, a <~ 1 
mb- and the px distribution of the 7r°'s. We have proceeded as follows: 
In pp collisions at y/s = 200 GeV, the shape of the pr distribution of 7r° 's can be 
described as (1+Pt/po)~™ with n = 9.99 andp = 1-219 GeV/c. The corresponding 
average pr is < pr >= 2po/{n — 3) = 0.349 GeV/c. In AuAu collisions we assume 
that the pr distribution of 7r°'s at each b is given by the same shape, (1 +pt/po)~ i \ 
with n — 9.99 and changing the scale p into po{b) =< pr >b (n — 3)/2, where 
< Pt >b is the average value of pr measured experimentally at the corresponding 
centrality. For central AuAu collisions (n part = 350 , b = 2), < pt > e xp is equal to 
0.453 GeV/c resulting in a p of 1.583 GeV/c. 

We take n as fixed, since at large pt the shadowing vanishes and the ratio Raa 
is independent of pr- Also, the experimental value of n in dAu is the same as in pp 
within errors. 

We can thus compute the ratio Raa in the absence of final state interaction 
(column I, Table 1). In this case Raa increases with pr, there is shadowing are low 
Pt and we have a Cronin effect that does not decrease at large pr- 

To these values we apply the correction due to the suppression factor S % o : 
First, we neglect the second term -the gain- in eqs. (3-4) (colunm II, Table 1). 
We find a general suppression in Raa, independent of pt- Raa increases slightly 
with pr and agrees with experimental data for pt > 5 GeV/c. At lower pr the 
result is significantly lower than the data. 

Table 1. Values of R a °uAu(Pt) for 
the 10 % most central collisions 
AuAu collisions at mid-rapidities 
(\y*\ < 0.35). Column I is the result 
obtained with no final state inter- 
action. The other columns include 
final state interaction with several 
ansatzs for the px shift. 



Pt 


I 


II 


III 


IV 


V 


VI 


0.5 


0.38 


0.05 






0.34 


0.38 


2 


0.90 


0.13 


0.08 


0.11 


0.20 


0.31 


5 


1.48 


0.21 


0.18 


0.19 


0.23 


0.25 


7 


1.69 


0.24 


0.24 


0.24 


0.24 


0.24 


10 


1.84 


0.27 


0.34 


0.30 


0.25 


0.24 



Then we introduce the second term -the gain- in eqs. (3-4) in two ways: 
Assuming that the pr shift of the 7r°, due to its interaction with the medium, 
is constant (two cases : 5pr = 0.5 GeV/c and Spr = 1-5 GeV/c) (columns III and 
IV), we obtain a slight increase of Raa at large pr, rather insensitive to the value 
of the shift, consistent with data. The problem at small pr remains. 
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If we assume that Spx oc (pt~ < Pt >b) (columns V and VI), in such a way 
that the factor S^o is 1 at px =< pt > as it should be -no suppression at low pt~, 
we obtain a slight decrease of Raa at large pt, consistent with data and an increase 
of Raa at low pt -due to the shift of large pt particles- resulting in an agreement 
with data. 




10 2 4 6 8 10 

p, (GeV/c) p T (GeV/c) 



Figure 1. Left: Values of RauAu(pt) f° r the 10 % most central collisions (lower line) and 
for peripheral (80-92 %) collisions (upper line) at mid-rapidities (\y*\ < 0.35), using the 
PT shift as Spt = (pt~ < PT >&) /20, (solid line). The dashed line is obtained using 
SpT ~ (pt~ < PT >&) ' /20 for pt < 7 GeV/c and pt constant for px > 7 GeV/c. Right: 
Values of Rj,Au(Pt) f° r minimum bias collsions at mid- rapidities (\y*\ < 0.35), using the 
p T shift as Sp T = (p T - < p T > h ) 15 /20. The data are from PHENIxISl. 

We turn next to minimum-bias dAu collisions at central rapidity. Here < pt >= 
0.39 GeV/c. With the same value of n as above (n = 9.99), this corresponds to 
po = 1.346. Calculating the ratio dAu to pp we obtain the result of Fig. 1. At for- 
ward rapidity, RdAu decreases as y increases due two effects: 

The first effect is basically due to energy-momentum conservation. It has been 
known for a long time in hadron-nucleus collisions at SPS energies -low pt 
"triangle"-. Its extreme form occurs in the hadron fragmentation region, where the 
yield of secondaries in collisions off a heavy nucleus is smaller than the correspond- 
ing yield in hadron-proton. This phenomenon is known as nuclear attenuation. It 
turns out, that, at RHIC energies, this effect produces a decrease of RdAu of about 
30 % between rf = and rf = 3.2. 

The second effect is the increase of the shadowing corrections in dAu with increas- 
ing y* . This produces a decrease of RdAu(pr) between 77* = and 77* = 3.2 of 
about 30 % for pions produced in minimum bias collisions. Therefore, we expect 
a suppression factor of about 1.7 between RdAuiPr) at 77* = and at rj* = 3.2, 
practically independent of pt ■ This is consistent with the BRAHMS results ^3. 
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Finally, there are some importants remarks to be done: 

• We use jet absorption, which is equivalent to jet quenching. It points to the 
existence of a dense system. 

• Our final interaction takes place at very short times -with a medium of 
partonic or prehadronic nature-. This is not a rescattering of fully formed 
hadrons. 

• We have studied the different possibilities of Spr- This is equivalent to a 
test on whether the mean energy loss -pr being E at y = 0- is constant 
with energy -which is the BDMPS ^ result for large enough energies-. 

• Our results for pt > 5 GeV/c depend little on the form of Spr- The in- 
crease of the ratio with increasing pt for a constant 5pr coincides with jet 
quenching results when reasonable kinematical constraints are imposed . 

• In perturbative QCD, Raa should tend to unity at large pr- However, 
this may occur at much larger values of px than the present ones. The 
BDMPS energy loss result being energy independent applies for parton 
energies larger than those available at RHIC. 

• We use Cronin effect in an effective way. Our Cronin does not die at large 
Pt, and it does not increase with y either. 

• For the region in which x is small enough, we use shadowing. It would 
lead to saturation at very small x. So our shadowing provides a transition 
between saturation at very small x and no shadowing at large x ~ 0.1. 
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